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EGO LAUNCH WINDOW STUDY 


% 


I. INTRODUCTION 

The purpose of this study is to determine whether a launch time 
exists, such that the practical restraints, imposed by the spacecraft 
and experiments, can be satisfied. This report deals primarily with 
the spacecraft imposed restraints. The experiment imposed restraints 
will be considered in a later report. 


II NOMENCLATURE 


* 


a = Semimajor axis of the elliptical orbit ~ kilometers 

e = Eccentricity of the orbit 

i = Inclination = The angle between the orbit plane and the 
equatorial plane ~ degrees 

O = Longitude of ascending node — the angular distance from 
•the vernal equinox measured eastward in the equatorial - 
plane to the point of intersection of the orbit plane where 
the satellite crosses from south to north ~ degrees 

Line of Apsides - The line connecting perigee and apogee 

u) = argument of perigee = the angular distance measured in the 

orbital plane from the line of nodes to the line of apsides 
~ degrees 

h = Perigee height = this is defined as the distance from the surface 
^ of the reference sphere (of radius equal to the Earth's equation- 
al radius) to the spacecraft measured al ?rg the radius vector when 
the spacecraft is nearest to the reference sphere ~ kilometers 

$ = Geocentric latitude = this is defined as the angle at the center 

of the Earth between the radius through a given point and the 
equatorial plane ~ degrees 

6 = Terrestrial longitude = the angular distance from the Greenwich 

meridian, measured eastward along the equator to the meridian 
plane ~ degrees 

h = Geocentric height = this is defined as the distance from the 
reference sphere (of radius equal to the Earth's equatorial 
radius) to the spacecraft measured along the radius vector 

v = The speed of the vehicle ~ ft/sec 
s 

A c = Azimuth = the compass heading measured clockwise from north ~ degrees 

E = Elevation Angle = the angle between the velocity vector and tne 
plane normal to the radius vector passing through the vehicle 
- degrees 
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III. ASSUMPTIONS • 


A. Initial Conditions 

The nominal geocentric polar coordinates (See Section II for definition) 
at injection with their ± 3 a tolerance are as follows: 

$ - -13*927 ± 0.284 degrees 

8 = 124.372 ± 0.492 degrees 

h = 282878.0 ± 14636.0 meters 

v 8 = 10541.0 ± 14.0 meters per second 

A c = 61.932 ± 0.323 degrees 

E = 1.29126 ± O.O 98 degrees 

B. Launch Window Restraints 

1. The perigee height will stay above its initial value for a one 
(l) year period. This eliminates the possibility of the satellite losing 
its energy, due to aerodynamic drag, and plunging into the Earth. 

2. The maximum eclipse time per orbit shall not exceed two (2) hours 
for a one (l) year period. For eclipse times of a longer duration the 
temperature of the solar array would go below an allowable limit. 

3. At least 262 minutes will elapse between injection and the first 
eclipse. There will be a drain on the battery power supply from lift-off 
to acquisition of attitude. In order for the batteries to have sufficient 
energy to sustain the first eclipse, it will be necessary for the spacecraft 
to spend at least 262 minutes in sunlight between injection and the first 
eclipse. 

4. The injection point will not lie in darkness. The attitude control 
system is not designed to acquire the spacecraft in darkness. 

* Mr. R. E. Russey of the 0G0 Project Office furnished these initial conditions 
and restraints. 
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C. Mathematical Model 


1. Perturbing forces included in the model — The perturbing 
forces due to the Sun, the Moon, and the gravitation field of 
the Earth (2nd, 3 r d and 4th harmonic) were included in the 
analysis for the construction of the "Launch Window" map. 

(Figure l). 

2. Perturbing forces not included in the model -- The perturb- 
ing forces due to the aerodynamic drag and the radiation pressure 
were not included in the analysis for the construction of the 
"Launch Window" map. (Figure l). 

Figures 2 through 8 show that radiation pressure and 
aerodynamic drag have a negligible effect on the EGO orbit 
for a launch time of July 15, 1963 at k hour UT. Each figure 
presents three ( 3 ) curves, two of which were determined with 
the Interplanetary Trajectory Program, and the third with the 
Launch Window Program. The difference between the two curves 
generated with the Interplanetary Trajectory Program is that 
one includes the effects of aerodynamic drag and radiation 
pressure and the other does not . The third curve, the one 
generated with the Launch Window Program is included on the 
figures to show the agreement between the programs. 

The aerodynamic drag coefficient is defined in the conventional 
manner 



i P u 2 a 


where U is the flight speed of the spacecraft relative to the atmosphere 
of density p and A is the cross-sectional area of the body normal to the 
flight direction. A realistic value of C^A = 170 square feet was assumed. 
The value of = 2.5 -*-0.5 I s recommended by Sehamberg (Reference l). 

The cross sectional area A for EGO varies between 10 square feet and 90 
square feet. 

The acceleration due to solar radiation is (1/c) K (A/M), where I 

is the solar energy flux, C is the velocity of light and K is a constant 

(0 = K = 2) whose value depends on the reflecting characteristics of the 

surface. (Reference 2). For a satellite which is transparent to the 

solar flux I, K = 0. For specular reflection from a flat plate, K = 2. 

The cross-sectional area A which is normal to the vector from the Sun 

to the satellite varies from 10 square feet to 90 square feet and the 

mass of the spacecraft is about 1000 pounds. The solar flux I = 2.00 

cal/cm 2 - minute, K = 1 and (A/M) = 0.17409 square centimeters per gram 

was assumed for this analysis or equivalently a perturbing acceleration 

of about 8.1 x 10~ 6 77^ o 

sec . 

Figure 9 presents another comparison between the Interplanetary 
Trajectory Program and the Launch Window Program for a launch date of 
July 15 , 1963 , at 0.4 hours OT. For this launch date the perigee is lower 
and aerodynamic drag should have a larger effect. This comparison lends 
further credence to the assumption of zero aerodynamic drag and zero 
radiation pressure used in the Launch Window Program. 
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IV. METHOD OF ANALYSIS 

Two IBM 7090 computer programs were usee to generate the results 
of this report. A discussion of the two programs follows: 

A. The Launch Window Program 

This program was used to establish that part of Figure 1 
corresponding to the restraints: i.e., perigee shall not go below its 
initial value for a one (l) year period and the maximum eclipse time 
per orbit shall not exceed two (2) hours for a one (l) year period. 

The input quantities to this program are the orbit elements 
a , e , i , o and <ju at the initial time t . The program calculates 
a, e, i, Q and a' as a function of time 7, influenced by the gravitational 
forces due to the Moon, the Sun and the Earth. The time per orbit that 
the satellite spends in the Earth’s shadow is also in the program’s 
output. In other words, the program calculates the shape of the orbit 
(a, e) and the orientation of the orbit in space (i, p and m) as a function 
of time t, but it does not calculate the position the spacecraft has in 
the orbit (true anomaly, etc.) 

With regard to perigee, the perigee altitude is calculated for 
a time when the spacecraft is not in general at perigee. This may be done 
because perigee height changes gradually over one period of the satellite's 
orbit. 

The time that the satellite spends in shadow is based on a set of 
elements (a, e, i, Q, and <») and the Sun’s position vector which are for 
a time r. In general the spacecraft Is at some arbitrary point at time 
t and not necessarily in shadow at that time. It is assumed that the 
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orbital elements and the Sun's position vary so slowly over one period 
of the satellite's orbit that the shadow time per orbit can be calculated 
by using tne aforementioned elements, the Sun's positions and Keplei *s 

law. 

The gravitational effects of tne Moon are computed using Halphen's 
method as modified by F. Musen to facilitate the numerical computation. 
(Reference 3). This method permits the numerical integration of the long 
range lunar effects over an interval of many years. The short period terms 
are removed from the disturbing function by averaging it over the mean 
anomalies of the Moon's orbit and the satellite's orbit. Comparison 
between the curves calculated with this program and the Interplanetary 
Trajectory Program, (Figures 2 , 5, 7 and 9) show that this program 
averages out the periodic Moon terms (with a period of 26 days), whereas 
the Interplanetary Trajectory Program includes the short periodic Moon 
terms. 

The gravitational effects of the Sun are computed by means of a series 
expansion of the disturbing function. (Reference 4 ). The short period terms 
(of period equal to the satellite's period) are removed by averaging the 
disturbing function over the mean anamoly of the satellite’s orbit. The 
periodic terms (the terms with a period of l82f days or the time required 
for the Earth to revolve half way around the Sun) are included in the 
calculations by this program. This is evidenced by Figures 2 , 5, 7, 9, 

10 and 12. 

The Earth's perturbing gravitational potential is represented by the 
second, third and fourth zonal harmonics (page of Reference 5 )* 
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The effects of periodic terms are not included since for the pur pc 'e of 
this work only the secular effects are needed. The shadow computation 
was developed by A. J. Smith of the Theoretical Division. A. J. Smith 
also did the programming of the basic orbital computations including 
Halphen's method. R. Devaney and H. Montgomery of the Theoretical 
Division adapted the basic program to OGO problems. 

The advantage of this program is its speed. It takes about 4^ 
minutes of computer time to compute a one year orbit using a five ( 5 ) 
day integration interval. 

B. The Interplanetary Trajectory Program, (Reference 6) 

This program was used to establish that part of the "Launch Window" 
map (Figure l) which corresponds to the restraint: at least 262 minutes 
shall elapse between injection and the first eclipse, and the injection 
point will not be in darkness. 

This program was also compared with the shadow calculation of the 
Launch Window Program for mutual verification. 

The input quantities to this program are the geocentric polar 

coordinates, 0 , 0, h, v , A , E, at the initial time. (Options exist 

s c 

to input the initial conditions in other coordinate systems). 

A variety of outputs are generated but only those outputs of 
interest to OGO will be mentioned. The program calculates the osculat- 
ing elements a, e, i, rlJ , q, and the mean anomaly M corresponding to any 
subsequent time t, influenced by effects of additional bodies, the aid, 

3rd and 4th zonal harmonics (these are functions of the satellite's 
declination) and the aid tesseral harmonics (these are functions of the 


right ascension of the satellite with respect to Greenwich) of the Earth, 
aerodynamic drag, triaxial lunar potential, thrust, and radiation pressure. 

The times per orbit that the satellite spends in the refracted Earth's 
shadow (umbra and penumbra) are also included in the program's output. 

This is a versatile program but requires correspondingly greater 
amounts of computer time. It takes about 1 to 1^ hours of computer time 
to compute a one year EGO type orbit, 

RESULTS AND DISCUSSION 

A. Based on Nominal Injection Conditions 

Figure 1 presents a map of all of the possible launch days and 
injection hours from May 30, 1963 to July 3 , 1964 . Only those launch 
times corresponding to the areas labeled "all restraints satisfied" will 
satisfy the four restraints innumerated in Section III. Those areas labeled 
"eclipse exceeds 2 hours per orbit T days after launch" satisfy restraints 
1, 3 and 4 for a whole year but satisfy restraint 2 only for T days. For 
example, August l8, 1963 at 0 hours U.T., satisfies all restraints for 
330 days as labeled on the map. All other areas on the map are excluded. 

Figure 1 shows that July 15 , 19^3 ab 4 hours U.T. , is a completely 
acceptable launch date. The data presented in figures 2 thru 8 are for 
this l aunfh date. These figures show good agreement between the two computer 
programs and that aerodynamic drag and radiation pressure may be ignored. 

For this launch date. Figures 2 through 8 are as follows: 


FIGURE 

ORDINATE 

ABSCISSA 

2 

Perigee height, 

Days after launch 

3 

Eclipse time 

Days after launch 

4 

Argument of perigee, w 

Days after launch 

5 

Eccentricity, e 

Days after launch 

6 

Semi-major axis, a 

Days after launch 

7 

Inclination angle, i 

Days after launch 

8 

Longitude of the ascending 
node, n 

Days after launch 


Figure 2 shows that the perigee altitude goes from 280 km to 
3000 km in 360 days for this launch date. The periodic Moon term (of 
period of 26 days) shown in the two Interplanetary Trajectory Program 
curves (but averaged out in the Launch Window Program) of Figure 2 has 
a half amplitude of about 30 km. This should be considered when inter- 
preting the results of the Launch Window Program if perigee is allowed 
to go below its initial value where aerodynamic drag becomes appreciable. 

Figure 3 presents eclipse time per orbit versus days after launch 
as calculated by the Launch Window Program and the Interplanetary Tra- 
jectory Program. This figure shows that eclipses occur for about 30 
days duration with a maximum eclipse time of 1.8 hours per orbit and 
elipsesalso occur for 117 days duration with a maximum eclipse time of 
0.44 hours per orbit. 

Figure 4 shows that the argument of perigee ou changes from 329*5 
degrees to 360 degrees in 360 days. 
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Figure 5 shows that tne eccentricity e changes from 0,8915 "to 
about 0.8475 in 360 days for this launch date. 

Figure 6 presents semi -major axis versus days after launch for 
this launch date. 

Figure 7 shows that the inclination changes from about ‘31.07 
degrees to 41.0 degrees in 360 days. 

Figure 8 shows that the longitude of the ascending node Q changes 
from 141 degrees to 118 degrees in 36 O days. 

Figure 9 presents perigee altitude versus days after launch for 
a point on the "Launch Window" map which lies on the lower perigee 
boundary, i.e., a launch time of 0.4 hours UT on July 15, 1963 * This 
point was deliberately chosen on the boundary as a realistically worst 
case to show the effect of neglecting aerodynamic drag (the aerodynamic 
drag is higher at lower altitudes). 

Figures 10 and 11 present perigee height versus days after launch 
and shadow time per orbit versus days after launch for a marginally 
acceptable launch date. Figure 10 shows that the perigee restraint is 
met but Figure 11 shows that a two hour eclipse occurs 32y days after 
launch which violates the "2 hour eclipse restraint". 

In order to demonstrate the sensitivity of the perigee boundary, 
a point was considered which is not too far removed from the boundary. 

* Figure 1 shows .that July 15, 1963 at 0 hours UT is unacceptable, but 
is very near the perigee boundary. 
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Figure 12 shows that perigee altitude becomes unacceptable after 
only 78 days for this launch date. This demonstrates that the perigee 
boundary is very sensitive with respect to launch time. 

B. Injection Tolerance Study 

Table I presents the effects of the 3 a injection tolerances 
on the boundaries of the "Launch Window" map for a launch date of July 15, 
1963 . 

Columns 3 and 4 of Table I are labeled "lover perigee boundary" 
and "upper perigee boundary" respectively. The "low perigee boundary" is 
the earliest launch hour (UT) of July 15, 1963 for which the perigee re- 
straint enumerated in paragraph B of Section III is satisfied. The 
"upper perigee boundary" is the latest launch hour for which the afore- 
mentioned perigee restraint i6 satisfied. Columns 5 and 6 of Table I 
are labeled "lower shadow boundary" and "upper shadow boundary"; these 
columns give the earliest and latest times respectively, for which the 
2 hour eclipse restraint per paragraph B of Section III is satisfied. 

The last row of Table I defines the boundaries for the nominal injection 
conditions. A comparison of this row with the others show that the 3a 
tolerance affect the launch window boundary by at most l/4 hour. 

* 

CONFIDENCE IH RESULTS 

The confidence in this analysis is established by comparing the 
results of the Launch Window Program with the results of the Inter- 
planetary Trajectory Program and with experimental data. 


A. Comparison of Computer Programs witn Each Other 


The two programs are based on different theories and developed 
completely independently. Figures 2 through 9 show that the two programs 
are in good agreement. Further discussion is given under paragraph C 
of Section HI. 

B. Comparison of Computer Results with Exp eri mental. Data 

Figure 13 presents variation of perigee height for the Echo 
Satellite. 35 ae observed variation, indicated by the circles, is compared 
with the computed variation for a period of more than 580 days. This 
single computation carried out for 5Q0 days shows remarkable agreement 
for a numerical integration program over such an extended number of 
orbits. Hiis comparison shows that accurate orbital data are generated 
with these two computer programs. Therefore, there is a high degree 
of confidence in the data presented in this report. 
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TABLE I 

EFFECTS OF INJECTION TOLERANCES ON LAUNCH 'WINDCM for July 15, 1963 . 

Lower Perigee Upper Perigee Lower Shadow Upper Shadow 
Boundary Boundary Boundary Boundary 















